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We apply a method based on a BCS-type approach to the extended Hubbard model on a square 
lattice to deal with the Y Ba2Cuz07-5 family of superconductors under pressure. The parameters of 
the tight-binding band are taken from experiments, and the coupling strength U and V are estimated 
by the zero pressure phase diagram {Tc x rih)- This scheme yields the nontrivial dependence of the 
superconductor critical temperature Tc as a function of the hole concentration nn in the Cu02 plane. 
With the assumption that the pressure P modifies the potential V and the on plane hole content 
rih, we can distinguish the charge transfer and the intrinsic contribution to Tc{P). We show that 
the changes on Tc{P) for the Y Ba2Cuz07 optimally doped compound at low pressures are almost 
entirely due to the intrinsic term. 
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I. INTRODUCTION 

Pressure experiments have played an important hole 
in understanding high-Tc superconductors (HTSC). Sev.- 
eral review papers have been written on this subjectETH. 
In short, the pressure investigation is a helpful tool to 
display the potential effects of the chemical pressure, 
thereby providing useful information for synthesis of re- 
lated compounds and, more important to the present 
work, to understand which are the structural parameters 
that influence the microscopic superconducting mecha- 
nism responsible for the HTSC. 

It is well known that the pressure effects are twofold: 
(1) It is documented that the applied pressure Pincreases 
the hole concentration rih on the Cu02 planesci, with a 
maximum rate of charge transfer approximately given by 
drih/dP = 0.02GPa"^ ThisJs called the pressure in- 
duced charge transfer (PICT)q term. (2) Another pres- 
sure effect is the increase of Tc for the optimum doped 
compound {uop) above the zero pressure value, which 
indicates another mechanism independent of the PICT, 
known as the "intrinsic" contribution termQ. Recently, 
a very large intrinsic term of dTc/dP = 6.8KGPa~^ 
has been measuredtl, however, its origin is still a mat- 
ter of research. Usually these two effects, the PICT and 
the intrinsic term, have to be combined in order to ac- 
count for the experimental data on a given family of com- 
pounds. It is very convenient to separate both contribu- 
tion because the PICT can be measured independently 
by thermal measurements, like the Hall coefficient. Some 
earlier works have estimated the magnitude of both 
fects through a phenomenological pressure expansiono'EI 
applied to the Y Ba2Cu30j-s family of compounds. 

A method using a BCS-type mean field approach with 
an extended Hubbard Hamiltoniaa. was introduced to de- 
rive the phase diagram Tc x n/J9'Ll. This procedure was 
used to deal with the pressure effects and it was as- 
sumed that, besides the hole content n^, the pressure also 



changes the magnitude V of the attractive potential, and 
ultimately changes the zero temperature superconduct- 
ing gap A(0)ll3. Thus, the critical temperature Tc{nh,P) 
was obtained as an expansion in powers of P, for a deter- 
mined Uh value, and each contribution was singled out, 
which is very convenient to interpret the experimental 
results. This method was applied to the experimental 
results of Hg 1201 and Hg 1223, with a fixed set of pa- 
rameters, with excellent aspeement with the underdoped 
and overdoped compoundal3. In a similar approach, also 
using the BCS method-and the extended Hubbard Hamil- 
tonian, Angilella et alll used the isothermal compressibil- 
ity tensor to estimate the change of Uh with the pres- 
sure and, through the derived Tc x Uh, they obtained the 
variation of the attractive potential V with the pressure 
P. In their calculations they have also taken into ac- 
count pressure dependence for the lattice parameter ai, 
so that the hopping integrals ti depends on the pressure 
too {ti — ti{ai{P))). In this way, they could separate the 
charge transfer and the intrinsic contribution. Their cal- 
culations for the intrinsic term were ija-good agreement 
with the measurements of the Bi 22120. 

Calculations based on the extended Hubbard Hamil- 
tonian provide a possible microscopic interpretation for 
the intrinsic term since it is related with the supercon- 
ducting interaction V. Through such calculations, one 
may extract some clues to the behavior of the fundamen- 
tal microscopic mechanism under the structural changes 
provided by the applied pressure, and some hints for the 
mechanism itself. 

Recently, we have combined both approaches of Ref.i'i 
in order to study the measured Pii— j< P curves of the 
Tlo.5Pbo.5Sr2Cai-xYxCu207 seriesEJ. A characteristic 
of this method is that, if the phase diagram is experimen- 
tally known, TdP) can be calculated for the entire family 
of different n/i compounds, and the value of drih/dP cal- 
culated for this specific family. Our calculated Tc{nh,P) 
for the Tlii,^Pho,^Sr2Cai-xYxCu20j series agreed well 
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with- the experiments for the extended-s wave symme- 
tryi. 

In a recently letter, Chen et allia have calculated 
the pressure effects on the nearly optimum doped 
YBa2Cu30r and underdoped Y Ba2GuiO^, using a 
BCS-type method. Using values of duh/dP averaged 
from different results of several groups, they obtained 
very good agreement with the experimental data. How- 
ever, calculating dTc/dP {P w 0) for different com- 
pounds of the YBa2Cu30T-s family, they found almost 
the same values for different potential parameters, which 
led them to conclude that the main contribution to 
Tc{nh,P) was from the PICT term. To check whether 
the YBa2Cu307 compound has an intrinsic contribution 
term or not we have performed calculations on the YBCO 
system using our method. 

In our method we use a BCS approach with the ex- 
tended Hubbard model to calculate the phase diagram 
of the Y Ba2C'u30T-s system, which is compared with 
the experimental data. Then we take into account the 
effects produced by the pressure through a change in 
the intersite potential: as first order approximation we 
write V{P) + AV{P). Here AV{P) = ciP, with ci 
being a constant independent of P that will be defined 
later. Since the structural changes for typical pressures 
are small, we neglect its effect on the hopping integrals. 
For the density of hole carriers we use the well know 
dependence of Uh with Fu. From this we estimate the 
critical temperature Tc{nhiP), V{P)) as an expansion in 
terms of the pressure P. Our results showed good agree- 
ment with the experimental data. 



II. THE PHASE DIAGRAM 

To develop the dynamics of the hole-type carriers in 
the Cu-0 planes, we -adapt a two dimension extended 
Hubbard HamiltonianEnlj'EJ in a square lattice of lattice 
parameter a 

H = - Ujcl^Cja + UY^ n,^n,i 

where tij is the nearest-neighbor and next-nearest- 
neighbor hopping integral between sites i and j ; U is the 
Coulomb on-site correlated repulsion and Vij is the at- 
tractive interaction between nearest-neighbor sites i and 
J- 

The transformation of Eq. (|^) to the momentum space 
leads to the appearance of the dispersion relation. In 
order to compare with YBCO system, we used a disper- 
sion relation within a tight-binding approximation, which 
may be compared with the one estimated from ARPES 
measurementala. Thus, 



ffk = —2ti{cos{kxa) + cos{kya)) + 4t2 cos{kxa) cos{kya) 
—2t3{cos{2kxa) + cos{2kya)) 
+At4 cos{2kxa) cos{2kya) 
—At5{cos{2kxa) cos(fcya) -I- 

cos{2kya) cos{kxa)) ~ n, (2) 

where it was considered identical hopping integrals along 
both directions in the Cu-0 planes for the nearest- 
neighbor {tx~ty~ti), a different one for the next-nearest- 
neighbor (^2), and so on. Here it was considered a hop- 
ping from the first to the fifth neighbor. Equatiop-|(0) is 
in agreement with the one given by Schabel et all3; /i is 
the chemical potential which controls the hole concentra- 
tion. 

Like the low temperature superconductors, the HTSC 
exhibit an energy "gap" (order parameter) in the su- 
perconducting phase which separates the paired states 
from the single-particJc states. Using a BCS-type mean- 
field approximationtZl to develop Eq.(^ in the momen- 
tum space, one obtaiiia.the self-consistent gap equation, 
at finite temperaturesEZi 

A.--E^-~tanh^, (3) 

k' 

with 

= ^el + Al (4) 

which contains the interaction potential I4;k' that comes 
from the, transformation to the momentum space of 
Eq.(|ll)Be 

Vkk' = U + 2V cos{kxa) cos{k'^a) + 2V cos{kya) cos(fc^a), 

(5) 

Phase Diagram 
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FIG. 1. Phase diagram Tc x Uh for the d wave and 
extended-s wave symmetries, with rttie experimental data of 
the Y Ba2Cu30r-s taken from Ref.E3 (open squares). 
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The substitution of Eq. (||) into Eq. (||) leads to appear- 
ance of a gap with two distinct symmetriesa: 

Ak(M,T) = A"'^"(^,T)[cos(fc,a) ±cos(fcya)], (6) 

where the plus sign is for extended-s wave and the 
nus sign, for d wave symmetry. In accordance with Re: 
one observes that the d wave part of the gap do not de- 
pend on the couphng constant U, depending only on V. 
The extended-s sy-cametry depends on both U and V. 
As it is well knownEJ'EEl at the critical temperature there 
is no symmetry mixture, and the gap might be in the 
extended-s or d wave state when T Tc- 

Using the same BCS-type mean-field approximation 
used in ihe gap equation one obtains the hole-content 
equationEd 



k 



1 



Ek ^ , £^k 

-— tanh 



(7) 



where < < 1. This equation, together with the gap 
equatios and the dispersion relation, are solved numeri- 
cally self-consistently, in the limit of T — > in order to 
obtain the phase diagrams for both gap symmetry. 

To compare with the experimental results for the Y123 
system, we perform our calculations with the hopping in- 
tegrals: fi=0.32eV as the nearest-neighbor value; for t2 
it was adopted the ratio ^2/^1=0. 50 for the d-wave sym- 
metry and ^2/^1 =0.57 for the extended-s; the remaining 
band parameters followed: t^/ti—O.W, t4/ti=0.11, and 
^5/^1=0. 031. AU these valueSj-aj-e close to the ARPES 
measurements on Y123 systemEEl. 

Using these band parameters we first estimated the 
coupling constant V, using the d-wave gap symmetcjf, 
reproducing the experimental phase diagram of Ref.c3. 
The best value obtained was y=-0.1152eV. For the case 
of the extended-s gap symmetry, the value C/=0.0096eV 
was obtained, using the same coupling constant V of 
the d-wave. The chemical potential /x is calculated self- 
consistently. Fig. 1 shows the results of our numerical caL. 
culations together with the experimental data of Ref.EHl 
for the Y Ba2Cu-iO'!-s series. 



III. THE METHOD 

To calculate the variations of Tc for a given compound 
with a certain value n^i, and under pressure P, we may 
use an expansion of Tc(n;i, P) in powers of P. Therefore, 



dT^ 
dP 



P 



1 rd^Tr 



where 
d'^Tr 



dP' 



2! V dP'^ 



'^^ dV ^"^ duh 



(8) 



T,K(P),\/(P)). (9) 



Here, ci = (dV/dP) and C2 = (drih/dP) are constant pa- 
rameters, which are determined fitting the experimental 
data for a given system. Tc{nh, 0) is the critical tempera- 
ture for P=0. Equation (H) can be written in a compact 
form, as 



T,{nh,P) = Y.' 



(10) 



with 



Cl 



dV 



C2 



dm 



The first coefficient (z 
jiven by 



r,K(p),y(p)). (11) 



1) for the above expression is 



ai 



Cl 



W 



duh 



(12) 



It is important to stress that the first term in ai is the 
intrinsic contribution, and the second one is the PICT 
contribution. Therefore, both contribution can be sin- 
gled out. Restricting ourselves to small changes V, we 
may approximate 



dTc 
dV 



AV 



(13) 



Phase Diagram 
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FIG. 2. Phase diagram Tc x rih for the d-wave showing 
the effect of a change on the coupling constant V, together 
witlj-the experimental data of the Y Ba2Cu^07-s taken from 
RefE3 (open squares). 
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where the horizontal bar denotes a "mean" over the 
values of Tc obtained direct from the phase diagram as a 
function oiV, as it is shown in Fig. 2. This "mean" is re- 
alized for each value of n/j. To obtain dTc/ drih we may ei- 
ther use an experimental value, or use directly the curves 
of Fig. 2, as done here. One can also use a phenomcnploft- 
ical universal parabolic relation between Tc and n;iE3'll3. 
For the coefficient a2 we have 



YBa^CUgO^. 



100 



T.(K) 



/ d d \ f dTc dTc\ 



80 - 



■ n„p = 0.165 
' n, =0.15 



Thus, a2 becomes 



70 



10 

p 



15 



20 



a2 



d'Tc 

dnl 



(15) 



The other coefficients of the expansion can be obtained 
using the same procedure. 



IV. THE PRESSURE EXPERIMENTAL DATA 



First of all, one observes that, at low pressures only 
the linear term (ai) of the expansion Tc{nh,P) comes 
into play. Therefore, we can approximate ai as the slope 
of the initial points of the Tc x P curve, for a given 
compound. Starting with nh=nop, we determine ci us- 
ing the estimated ai in Eg. (|l^) , as long as at Uop the 
charge transfer term vanishes. To determine C2 we ob- 
tain ai from the Tc x P curve of an Uh 7^ nop compound, 
and use again Eq.([l^). Once these two constants are 
determined, the coefficients for any rih value can be 
calculated (Eq.(|ll)). 

For the Y Ba2Cu^OT-s system we only have the non 
zero pressure data for the n^=rio^pp.l65 {5 ~ 0) and 
nh=0-l5 {6 w 0.15)|-C|(MHpoundsOo. Thus, from the 
TcXP curves of Ref.Bbl we estimate ai = 0.75K/GPa 
and 1.4 K/GPa for riop and ri;i=0.15, respectively. The 

phase diagram parameters (-^y-) taken from Fig. 2 and 
used in the calculations were lQ't>2K/eV and lAlbK/eV 
for the d-wave and llMK/eV and im%K/eV for the 
extended-s wave for Uop and n/i=0.165 compounds, re- 
spectively. Therefore, the resulting constants for the 
d-wave were ci = 4.54 x lQ~^eV/GPa which furnishes 
an intrinsic contribution of 0.8K/GPa, and C2 = 3.15 x 
10~^GPa~^ which is the charge transfer term (c2 — ^p-)- 
For the extended-s we obtained Ci = 6.4 x 10~'^eV/GPa 
and C2 = 2.2 x 10~^ G Pa~^ . It is important to men- 
tion that our values for dT^/dP (the intrinsic contribu- 
tion; see Eq.(p^) were in the|-same order of magnitude 
of Neumeier and ZimmermannO, despite the difference in 
the approaches. 



FIG. 3. Numerical results: the solid and dashed lines 
are the d-wave results, while the long-dashed lines are the 
extended-s wave results. The open squares and filled squaries 
are expsfjpiental data of Y Ba2Cu307-s taken from RefB3 
and Ref£j, respectively. 

On Fig. 3 we compare the Tc{nh, P) expansion with the 
experimental data. We observe a very good agreement 
for both compounds, specially for Uop, using the d-wave 
gap symmetry, where the data can be fitted up to 20 
GPa. The increase of riop with the pressure is attributed 
to the intrinsic contribution in our theory, and the maxi- 
mum Tc at P ~ 7GPa and subsequent decrease is due to 
the negative second order term of the expansion, which 
reflects the effect of the PICT. The extended-s wave re- 
sults did not show good results indicating a preference 
for the d-wave pairing mechanism for the Y-123 system. 



V. CONCLUSIONS 



We conclude this letter emphasizing that our method 
based on the BCS approach with the extended Hubbard 
model, with both the attractive potential V and the den- 
sity of carriers Uh depending on the pressure P, predicted 
an intrinsic term and aj-PICT contribution similar to 
those predicted earlier BlO. In our method we propose 
a direct pressure dependence for the attMctive potential 
V, which is different from other worksa that use V as 
an adjustable parameter in order to fit the known exper- 
imental dependence of Tc on P. Our results for the d- 
wave are in excellent agreement with the Y-123 pressure 
data and, the intrinsic term being the most important 
contribution to the low pressure data on the optimum 
compound. Concerning the underdoped compound the 
intrinsic and PICT contribution are both present and 
must, aajpreviously known by several experiments and 
theorieffliJ, be taken into account. Such results support 
the d-wave symmetry in these compounds. 

We thank CNPq, Capes and FAPERJ for finacial sup- 
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